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Abstract 

We investigate the optical properties of edge-functionalized graphene nanosystems, focus- 
ing on the formation of junctions and charge transfer excitons. We consider a class of graphene 
structures which combine the main electronic features of graphene with the wide tunability of 
large polycyclic aromatic hydrocarbons. By investigating prototypical ribbon-like systems, we 
show that, upon convenient choice of functional groups, low energy excitations with remark- 
able charge transfer character and large oscillator strength are obtained. These properties can 
be further modulated through an appropriate width variation, thus spanning a wide range in 
the low-energy region of the UV-Vis spectra. Our results are relevant in view of designing 
all-graphene optoelectronic nanodevices, which take advantage of the versatility of molecular 
functionalization, together with the stabiUty and the electronic properties of graphene nanos- 
tructures. 

Keywords: carbon nanostructures, configuration interaction, ZINDO, exciton, UV-Vis spec- 
trum 
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Research on the graphene chemistry has been stimulated in the last few years by an increasing 
demand for both large scale production of graphene samples and atomic control on its nanostruc- 
tures. The electronic features of graphene, related to its ;r-electron network, make this novel 
material analogous to large polycyclic aromatic hydrocarbons (PAHs)^'-^ and allow to exploit a 
wide range of well-known functionalizations, reactions and preparation techniques, also useful for 
technological applications which need the blending of graphene with organic and inorganic com- 
pounds. One of the main advantages of these chemical approaches is related to the possibility 
to obtain graphene nanostructures, tailoring their size, shape and edge termination on the atomic 
scale. "^'^ The control of these parameters can be exploited to tune the electronic ^"^ and optical 
properties of graphene nanostructures. 
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Figure 1: Scheme of the energy offset between edge covalently functionalized graphene 
nanoflakes, with respect to a hydrogenated flake (H-flake) taken as a reference. Edge functional 
groups with remarkable electron-donating (-withdrawing) character, such as -NH2 (F) induce an 
almost rigid upshift (downshift) of the gap region. 

In this direction, we recently demonstrated^^ that edge covalent functionalization represents a 
viable method to tune the ionization potential (IP) of graphene nanoflakes (GNFs), preserving the 
overall characteristics of the ;r-conjugation. Furthermore, as illustrated in Figure 1, a correspond- 
ing shift of the electron affinity (EA) moves the energy gap region of the GNF in an almost rigid 
fashion according to the edge termination. With respect to a reference hydrogenated flake (see 
Figure 1), an upshift comes from electron-donating and a downshift from electron-withdrawing 
functional groups. The effectiveness of this mechanism allows us to form all-graphene type II 
nanojunctions (GNJs) which can be built such that occupied and virtual molecular orbitals (MOs) 
close to the gap region localize on opposite sides. This opens the possibility of obtaining charge 
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transfer excitons, which are relevant for many optoelectronic applications. 

Here we study the optical properties of such edge-functionalized GNJs, with particular atten- 
tion to the conditions for the formation of charge transfer excitons. We start characterizing the 
UV-Vis spectra of hydrogenated and functionalized GNFs, as obtained by considering atomic flu- 
orine (F), methyl-ketone (-COCH3) and amino (-NH2) groups. We then focus on few prototypical 
nanoj unctions, obtained by varying in turn edge functionalization and junction width. We show 
that upon appropriate modulation of these parameters, GNJs can be engineered in order to get large 
oscillator strength and charge transfer character for the first excitation. Moreover, following the 
energy gap modulation with respect to the width, the energy of the first peak in the UV-Vis spectra 
can be tuned in order to span a wide range (~ 1 .5 eV) in the low energy region. 

The UV-Vis spectra are computed within the framework of the semi-empirical Hartree-Fock 
(HF) based method ZINDO, through the configuration interaction (CI) procedure, including 
single excitations only (ZINDO/CIS). This method is known to provide reliable results for aro- 
matic molecules. ^^'^^ All calculations are performed starting from optimized geometries (0.4 kcal • 
mol^VA force threshold) and mean-field ground state properties calculated with AMI. 
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Figure 2: UV-Vis spectra of functionalized graphene flakes (empty solid curves) and of the ref- 
erence hydrogenated flake (shaded gray area). The curves represent oscillator strength (OS), nor- 
malized with respect to the number of C rings, vs. excitation energy and are convoluted with a 
Lorentzian function, with a broadening of 0.05 eV. The structure of the hydrogenated (H-) flake 
and the reference axes are also shown. 

We first discuss the optical properties of homogeneously functionalized graphene flakes in 
comparison with those of their hydrogenated counterpart (the structure is shown in Figure 2). 
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These systems are characterized by armchair-shaped edges and width parameter N =S, where 
indicates the number of dimer lines along the zigzag direction, borrowing the standard notation 
used for quasi- ID armchair graphene nanoribbons (AGNRs).^^ We have already shown that the 
class of graphene flakes considered here follows the width-modulated behavior of the energy gap 
(Eg equal to the difference between IP and EA) typical of AGNRs, also upon edge functional- 
ization. In particular, while the band gap of AGNRs scales inversely with respect to the ribbon 
width w,^^'^^ three ~ 1/w scaling laws have been identified depending on the parameter, ex- 
pressed as A'^ = 3;? + m (/? is a positive integer and m = 0, 1 , 2). The smallest gap pertains to the 
N — 3p + 2 family. The flakes considered here are about 3 nm long (jc direction) and their borders 
along the width (y direction) are shaped in order to minimize their influence on molecular orbitals 
(MOs) and to prevent the presence of localized magnetic states. The functionalizations treated 
here comprise both -NH2 groups, which induce an upshift of the gap region with respect to the 
hydrogenated counterpart, and -COCH3 groups and F atoms, which on the contrary downshift the 
gap region, according to Figure 1. The synthesis of fluorine substituted graphene molecules ^'^^ 
as well as functionalized with -NHj^ and ketone-based^'^ groups is well known and has been 
experimentally demonstrated. COCH3- and NHj -functionalized flakes are here obtained by re- 
placing each second H-terminating atom along both edges with a covalently bonded group, while 
atomic fluorine substitutes each H edge atom in the case of F-flake (see upper panel of Figure 1). 
Both functionalized flakes and junctions display distortions due to steric effects induced by edge 
terminating groups. 

By inspecting the UV-Vis spectra of the flakes shown in Figure 2, we notice that the optical 
properties of the systems are not significantly changed upon functionalization, at least in the low 
energy window. Only a slight red shift is noticed, on account of the reduced energy gap with respect 
to the hydrogenated flake: this effect is maximized in case of NH2- and F- functionalization, which 
show almost equal electrical and optical gaps. The first excitation, labelled as El (8) in Figure 2, 
is optically active in all cases and dominated by a HOMO LUMO transition. The second 
peak in the spectra [E3(8)] is generated by the third excitation and is mainly due to HOMO-1 
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LUMO and HOMO — )■ LUMO+1 transitions. Although these transitions are allowed by symmetry, 
the overlap between the involved orbitals is small and thus the oscillator strength for E3(8) is 
reduced of about one order of magnitude with respect to El (8). Between the two we find another 
excitation [E2(8)], which is very weak but not forbidden by symmetry, thanks to a small but non 
vanishing contribution of HOMO — )■ LUMO+1 transition. Further details on the composition of the 
excitations and on the symmetry of the involved MOs are reported in the Supporting Information. 




0.0 0.5 1 1.5 2 2.5 3.0 3.5 

Energy [eV] 



Figure 3: (a) UV-Vis spectrum of A'^ = 8 H//COCH3 nanojunction. The curve represents the oscil- 
lator strength (OS), normalized with respect to the number of C rings, vs. excitation energy and 
is convoluted with a Lorentzian function, of broadening 0.05 eV. The structure of the junction is 
also shown; (b) Electron (e) and hole (h) probability density for E1(J8) excitation of H//COCH3 
junction, after Eq. (3) and Eq. (2). 

We next consider the optical properties of the nanojunction obtained by interfacing H-termination 
and COCHj-functionalization: this system has the same width of the GNFs discussed above (i.e. 
N = S) while its length (~ 4 nm) is tailored in such a way that each side of the junction cor- 
rectly accounts for the properties of the corresponding homogeneously functionahzed system. By 
inspecting the UV-Vis spectrum of this junction, shown in Figure 3(a), we notice that the shape 
resembles that of isolated GNFs. However, the character of the first excitation, labelled as E1(J8) 
in Figure 3(a), changes dramatically. In fact, HOMO and LUMO of the H//COCH3 GNJ, which 
mainly contribute to E1(J8), are localized on opposite sides of the junction (see Fig. S2 in the 
Supporting Information). A charge transfer character is thus expected for this excitation. 

In order to better evaluate the charge transfer character of the excitations, we compute the 
spatial distribution of their hole (h) and electron (e) components. ^^'^^ Taking advantage of the 
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standard LCAO expansion of MOs^" adopted here. 



^i{r)^Y^ajiXj{r), (1) 
j 

where Xji''') Slater-type orbital basis functions and ap are the projection coefficients, we ex- 
press the probability density of h and e for the P'^ excited state as: 

P^(r) = Il4^-0a(r)P (2) 

pi(r) = E|cl,^-</)Mr)|2, (3) 

where c'^^ are CI coefficients weighting the contribution to the I''^ excitation of each transition 
a — > j8 from occupied (0a) to virtual (<j)p) MOs. The spatial localization of each quasiparticle, h 
and e, can then be quantified by computing: 

, Ly,zLx<oPh{x,y.z) 

= — —r X— (4) 

J __ I.y,zI.x>()Pe{x,y.z) 
Y.x,y,zPe\X,y.z) 

where the reference axes are centered at the flake, as shown in Figure 3(a). In the case of a ho- 
mogeneous distribution, L/j/g amounts to 50% of the charge density, while it is equal to 100% for 
a fully locaUzed component. For comparison, we also estimate the spatial localization of frontier 
orbitals, HOMO and LUMO, which mainly contribute to the first optically active excitation (see 
Supporting Information, Table SII). The isosurfaces representing the hole and electron contribu- 
tions to E1(J8) are shown in Figure 3(b) and the corresponding values L/j and are reported in the 
third column of Table 1. We notice that, while the first excitation indeed presents a partial charge 
transfer character, the spatial localization of h and e is rather scarce. This feature, which is clearly 
inherited by frontier orbitals (see Table 1), is also accompanied by a large oscillator strength of 
E1(J8), as a result of the non negligible overlap of the orbitals contributing to the excitation. 



Table 1: Spatial localization (per cent) of hole {Lh) and electron (Lg) components of the first 
optically active excitation of the H//COCH3 and NHj/ZF GNJs, belonging to the families = 
6, = 7 and = 8. The same for HOMO (Lhomo) and LUMO (Llumo)- 





H//COCH3 


NH2//F 


N 


6 7 8 


6 7 8 


Lh 

Le 


67 % 63 % 56 % 

68 % 75 % 60 % 


77 % 87 % 64 % 
87 % 74 % 65 % 


Lhomo 
Llumo 


80 % 73 % 58 % 
79 % 86 % 62 % 


92 % 96 % 66 % 
96 % 92 % 68 % 



In order to enhance the charge transfer character of the excitation, we exploit the potential of 
edge functionalization, which allows to maximize the energy offset at the interface by combining 
appropriate functional groups according to the scheme summarized in Figure 1. Moreover, we 
investigate the generality of this mechanism by varying the junction width over the three families, 
i.e. N = 3p, N = 3p+ I, and N = 3p + 2, which are known to present different electronic prop- 
erties also upon functionalization. In the following we consider prototypical junctions obtained 
by functionalizing one side with amino groups (-NHj) and the other side with atomic fluorine (F). 
The former have a very strong electron-donating character, as a consequence of the lone pair of 
N, which is shared in the bond with C edge atoms. The remarkable polarization of this covalent 
bond, with about -0.5 e on N and +0. 1 e on C - as estimated through MuUiken charge analysis 
- induces a spreading of frontier orbitals also onto N atoms. Noteworthy -NHj termination is pre- 
dicted to be quite stable and effective in modifying the electronic properties of quasi- ID graphene 
nanoribbons.-^^^-^^ On the other hand, F is strongly electronegative, on account of its open shell 
electronic configuration, typical of halogen atoms. Analogously to H, F atoms form g bonds at 
the edge: thanks to this similarity F-termination has already been proposed to tune the electronic 
properties of graphene nanostructures. 14,35,36 

We investigate the optical properties of NHj/ZF nanojunctions of width parameters N = 6, 
N = 1 and N =S. The optical spectra and hole and electron densities for the first optically active 
excitation are displayed in Figure 4; the spatial localizations are reported in the last three columns 
of Table 1. For completeness, we also include in Table 1 the spatial localization for the corre- 
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spending excitation in H//COCH3 junctions; their probability densities are shown Figure S5 of the 
Supporting Information. As can be seen from Figure 4, the UV-Vis spectrum of = 8 NH2//F 
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Figure 4: UV-Vis spectra of width modulated NH2//F graphene nanojunctions, belonging to fam- 
ilies N — 6, N — 1 and N = S. The curves represent oscillator strength (OS), normalized with 
respect to the number of C rings vs. excitation energy and are convoluted with a Lorentzian func- 
tion, with a broadening of 0.05 eV. Electron (e) and hole (h) probability densities are also shown 
for the first optically active excitations of each family. 

junction retains the basic features already discussed for the H//COCH3 functionalization, but is red 
shifted by 0.25 eV as a consequence of the larger energy offset at the junction interface, which 
reduces the effective optical gap. Consistently with the larger energy gap typical of the other two 
families, the first peak shows up instead at 2.09 and 2.55 eV for = 7 and N = 6, respectively. 
The comparison of the three GNJs shows that the width parameter can be engineered to tune the 
energy of the first excitation, which is always bright and intense, over a wide range of about 1 .5 eV. 
The effectiveness of the new NHj/ZF functionalization for the enhancement of the charge transfer 
character of the first optically active excitation is evident in the comparison of both probability den- 
sities (Figure 4) and spatial localizations of electron and hole (Table 1), the latter increasing up to 
20% with respect to the H//COCH3 case. Additional calculations performed on longer graphene 
nanojunctions (up to 5.5 nm) for the same NHj/ZF functionalization indicate equal or further in- 
creased spatial localization, confirming the validity and generality of our results. It is finally worth 
noting that in the N = S case the spatial localization of hole and electron is again very similar to 
that of the frontier orbitals, which mainly contribute to the first excitation. This is not the case for 
N = 6 and N = 7, where HOMO — > LUMO transitions represent only the 20% of the total weight. 
Nonetheless, additional transitions involving higher-energy localized MOs contribute to the first 
excitation in the N — 6 and N — 1 GNJs, thus leading to a prominent charge transfer character of 
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the first excitation. 

In summary, our results show that low energy excitations with remarkable charge transfer char- 
acter and large oscillator strength can be obtained for edge-functionalized graphene nanojunctions. 
Both prototypical H//COCH3 and NHj/ZF nanojunctions have been analyzed. A convenient choice 
of functional groups, which gives rise to opposite gap region shifts relative to the H-terminated 
flake, helps maximizing the energy offset between the two sides of the junction, thus enhancing 
the charge transfer character of the first optically active excitation. Furthermore, it is here demon- 
strated that a wide energy range (~ 1.5 eV) can be spanned in the optical spectra by means of an 
appropriate width modulation, retaining at the same time the charge transfer character of the exci- 
tation. The versatiUty of edge covalent functionalization, combined with the possibility of varying 
the width, makes these systems unique, opening promising avenues for designing all-graphene 
optoelectronic nanodevices. 
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